Quantum Size Effect in Conductivity of Multilayer Metal Films 
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Conductivity of quantized multilayer metal films is analyzed with an emphasis on scattering by 
rough interlayer interfaces. Three different types of quantum size effect (QSE) in conductivity are 
predicted. Two of these QSE are similar to those in films with scattering by rough walls. The third 
type of QSE is unique and is observed only for certain positions of the interface. The corresponding 
peaks in conductivity are very narrow and high with a finite cutoff which is due only to some other 
scattering mechanism or the smearing of the interface. There are two classes of these geometric 
resonances. Some of the resonance positions of the interface are universal and do not depend on the 
strength of the interface potential while the others are sensitive to this potential. This geometric 
QSE gradually disappears with an increase in the width of the interlayer potential barrier. 

PACS numbers: 73.40.-c, 73.63.Hs, 72.10.Fk, 73.50.-h 



I. INTRODUCTION 

Boundary scattering is essential for complete descrip- 
tion of nanosystems such as quantum wells, ultrathin 
films or wires, etc. Due to the large surface-to- volume ra- 
tio, boundaries are expected to play a much greater role 
in determining the overall properties in a nanostructure 
than in a bulk material. For example, recent scanning 
tunneling microscopy (STM) data have shown that elec- 
tron energy spectra can be more strongly correlated to 
the buried interfacial lattices than to the surface immedi- 
ately beneath the STM tip . These observations clearly 
indicate that a small lateral variation along the boundary 
can have a significant long range effect in a semi-ballistic 
electron system. Thus, a more realistic description of 
a nanoscale quantized system must go beyond the com- 
mon perfect geometric boundary and include boundary 
corrugations. Indeed, random surface roughness of a thin 
metal film can dominate incoherent scattering and relax- 
ation, and can lead to anomalous quantum size effect 
such as large oscillatory dependence of the in-plane con- 
ductivity on the film thickness Q{. 

The same must be true not only for the quantum well 
(film) walls but also for the interlayer interfaces in mul- 
tilayer films. It is well known that the roughness of the 
interlayer interfaces plays an important role in, for ex- 
ample, giant magnetoresistance (see review |3| and ref- 
erences therein). The purpose of this paper is to ana- 
lyze the effect of irregular corrugation of the interlayer 
interfaces on the lateral conductivity of quantized mul- 
tilayer films without magnetic effects. We will see that 
the interface scattering can result in unique features of 
the quantum size effect (QSE) which are strikingly dif- 
ferent from QSE with scattering by bulk or wall inhomo- 
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geneities. Orbital and spin magnetic effects of the type 
studied in Refs. |^ will be studied separately. 

In ultrathin films, the motion of electrons across the 
films can be quantized. QSE in metal films is studied ex- 
perimentally by measuring conductivity yilg] and suscep- 
tibility of the films or in spectroscopy Qand STM Q 
measurements (for earlier results, see references therein). 
As a result of QSE, the 3D electron spectrum e (p) splits 
into a set of minibands 6j (q) where q is the 2D momen- 
tum along the film (yz— plane). In the simplest case of 
a single-layer film approximated by a rectangular quan- 
tum well, the quantized values of the x— component of 
momentum are p Xj = nj/L (here and below % = 1). If in 
such quantized metal films the Fermi energy Ep is un- 
affected by the quantization, the Fermi surface reduces 
to a set of 2D curves e F j (q) that correspond to cross- 
sections of the 3D Fermi surface e (p) = Ep by a set of 
planes p Xj = irj/L, e Fj (q) = e F (p Xj ,<l)- 

This quantization of motion, which is determined by 
the film thickness L, leads to several types of QSE. First, 
any change of the film thickness L results in the change 
in size and number of the Fermi curves e F j (q). This 
thickness-driven change in number of the Fermi curves 
tpj (q) [or, what is the same, number of occupied mini- 
bands ej (q)] leads to a singularity in the density of states. 
These singularities are the most obvious manifestations 
of QSE. 

These singularities in the density of states, by itself, 
do not lead to any striking anomalies in the dependence 
of the lateral conductivity a of the film on the thickness 
L. The conductivity is more sensitive to electron scat- 
tering than to the density of states. However, the change 
in the number of the occupied minibands S can be ac- 
companied by a change in number of allowed scattering 
channels that correspond to the scattering-driven elec- 
tron transitions between minibands tj (q). The effect of 
this step-like change in the number of scattering chan- 
nels on the conductivity is much stronger than that of 
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the singularities in the density of states [9j- When all 
scattering-driven interband transitions are allowed, QSE 
manifests itself as a pronounced saw-like dependence of 
the conductivity on the film thickness. This type of QSE 
in quantized films has been predicted both for scattering 
by impurities and surface inhomogeneities [ToL Hll |. 

When the main scattering mechanism is the scattering 
by surface inhomogeneities, many of the interband tran- 
sitions can often be suppressed. This happens, for exam- 
ple, when the average size of the surface inhomogeneities 
R is much larger than the the thickness of the film and/or 
the particle wavelength Xp. Then the usual QSE, which 
is described above, disappears and is replaced by a differ- 
ent kind of the size effect Q. This anomalous QSE, which 
is somewhat reminiscent of the magnetic breakthrough, is 
completely decoupled from the singularities in the density 
of states and is associated solely with opening of inter- 
band scattering channels for gliding electrons at certain 
values of the film thickness, Lj ~ + 1/2) RXp/2. 

The main goal of this paper is to analyze QSE in mul- 
tilayer films with an emphasis on the scattering by the 
interface between the layers. We will see that, in addi- 
tion to the above two types of QSE, the multilayer films 
can exhibit a peculiar "geometric" QSE with very narrow 
high peaks in the lateral conductivity. Some of the posi- 
tions of these spikes in conductivity are universal; these 
spikes appear when the ratio of the thicknesses of the film 
layers is given by simple fractions. The position of the 
rest of the spikes depend on the strength of the interlayer 
interface. 

In the next Section, we briefly present the main equa- 
tions for the conductivity and introduce proper dimen- 
sionless variables. The results are presented in Section 
III. Section IV contains the summary and brief discussion 
of the results. Appendix contains auxiliary information 
on the energy spectrum of multilayer films of the type 
used in the calculations. 

II. CONDUCTIVITY 
A. Scattering by the interlayer interface 

For simplicity, we consider an ultrathin film of thick- 
ness L consisting of only two layers with the thickness of 
Li and Li- The interface between the layers is rough with 
random corrugation. The exact position of the interface, 
x = L\ + £ (y, z) , is described by the random function 
£ (y, z) with the zero average, (£) = 0. The random in- 
terface inhomogeneities £ (y, z) are best characterized by 
the correlation function Q (s) , 

C(s) = c(M) = «(si)£( Sl + s)> 

= a- x /^(soetsi+fOdsi, (i) 

where the vector s gives the 2D coordinates along the 
interface and A is the averaging area. Here, it is as- 
sumed that the correlation properties of the surface do 



not depend on direction. Two main characteristics of the 
surface correlation functions £ are the average amplitude 
(" height" ) and the correlation radius (" size" ) of surface 
inhomogeneities, t and R. 

To emphasize the scattering by inhomogeneities of the 
interlayer interface, we start from films with ideal outside 
walls that do not contribute to electron scattering. The 
combined effect of interface and wall inhomogeneities will 
be considered elsewhere. 

Mostly we are interested in the dependence of the lat- 
eral conductivity on the film thickness and have in mind 
the following experimental situation. The first layer of 
the film is grown on some (ideal) substrate. The surface 
is then roughened by adding inhomogeneous adsorbate 
or by some other means. The growth of the second layer 
starts from this rough interface and the conductivity is 
measured at different values of L 2 either in the process of 
growth or after the growth is completed. An advantage 
of such setup with a buried interface is that it allows one 
to measure the conductivity at various values of the film 
thickness with exactly the same random rough interface. 

In this setup, the thickness of the first layer, L\, should 
be considered as fixed, while the thickness of the sec- 
ond layer, hi, is variable. Below we are calculate the 
film conductivity a as a function of the film thickness 
L = L± + L2, c(L), assuming that L\ = const. The 
measurements of conductivity can be performed in sta- 
tionary conditions at different values of L2 or as a func- 
tion of time, in the process of film growth as in Ref. |12| . 

The second layer can be made of the same or different 
material as the first. If the material is different, then 
the electron potential energy between the layers differs 
by some AU. The the structure of the energy spectrum 
becomes a complicated function of AU making the be- 
havior of conductivity highly irregular |13| . 

Below we consider both layers to be made of the same 
material with the interface being the only disruption in 
the potential relief. Then the simplest model of the in- 
terface is the S~ functional potential barrier 

U = U 5(x-L 1 -£(y,z)). (2) 

This immediately introduces two new physical param- 
eters into the problem: the strength of the barrier Uq 
and its (average) position L±. In what follows, we study 
the dependence of the conductivity on these parameters. 
When necessary, instead of the (^-function we will study 
the corrugated interface with the finite width D. In 
experiment, the barrier can be a dislocation wall, twin 
boundary, or an oxide or dielectric layer (see, e.g., Ref. 
[T^ and references therein). 

The presence of the interface changes the spectrum. 
When calculating the changes in the spectrum, one can 
ignore small corrugation £ (y, z). The changes in spec- 
trum caused by the 8— type barrier (0) are discussed in 
Appendix. The random corrugation of the interface is 
responsible for the electron scattering and gives rise to 
the collision operator in the transport equation. 
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The scattering by the interface inhomogeneities leads 
to the transitions between the states ej (q) — > Cj (q'). 
Several ways of calculating the corrugation-driven tran- 
sition probabilities Wij (q, q') are described in Ref. fl3| . 
The simplest methods are either the direct perturbation 
approach ^| or the mapping transformation method 0] 
with both giving the same result in most of the parameter 
range. 

The corrugation-driven contribution 5U to the inter- 
face potential, Eq. © , with small corrugation £ is 



5U = -Uo£(y,z)5' (x-Lx). 



(3) 



The matrix element Vij (q, q') of this perturbation be- 
tween the states 6j (q) , tj (q') is 



=-U, 



o J e^-^i (s) % (x) 8' (x - L x ) *j (x) dxds 
Uot (q - q') [*< (Li) *5 (Li) + (L x ) (Lx)] , 



(4) 



where ^ (x) are the quantized wave functions for elec- 
tron motion across the film. Note that the derivatives 
1 J r ' (x) for films with a S— type barrier inside are discon- 
tinuous at the position of the barrier x = L\. Therefore, 
(L{) in Eq. (@J should be understood as % (Lx) = 
m(Lx + 0) + *' i (Lx-0)}/2. 

The corrugation-driven transition probability 
Wij (q, q') is given by the square of this matrix el- 
ement which should be averaged over the random 
inhomogeneities £: 

Wij (q,q') = (|Vy (q,q')| 2 ) ? = Etf< (|<fe-q$|) Gy, (5) 
Gij = [*< (Lx) (Lx) + V'i (Lx) (Lx)} 2 , (6) 

where Q (|qi— q^l) is the Fourier image of the correlation 
function of the interface inhomogeneities 0. The co- 
efficients Gik are calculated with the help of the wave 
functions presented in the Appendix. The explicit form 
of Gik is given in the next subsection. 

The transport equation is a set of equations for the 
electron distribution functions 7ij (q) in minibands and 
has the standard Boltzmann-Waldmann-Snidcr form 113: 



dru 
~dt 



2,Y j JWij[nj-ni]8(e itl 



" q, V) 2 ' 



The integration over dq' is done using the <5-function, 
<5 (e 4q - ejq') = m*jd (q' - q lj ) /q lj , where q lj (q) is the 
solution of the equation tj (q^ ) = e, (q) and the effec- 
tive masses m*j = qij/ (dej/dq) \ q= qiJ ■ As always in the 
transport theory, the angular integration is eliminated by 
using the angular harmonics. The current is given by the 
first harmonic of the distribution = Vi the equation 
for which involves only the zeroth and first harmonics 



W, 



(0,1) 



(q, q^) of W (q - q^ ) over the angle qq 
dUi (q)/dt=-^2 Vj (q^) /t^, 

3 

- L = fE[^ ) -^< 

W^=U^ :q q,U/, ; . 



(8) 



where we, to simplify the equations, assume that the ef- 
fective mass m*j does not depend on its indices, m = m*j . 

The solution of Eqs. (jSJ provides the 2D conductivity 
of the film: 



3h 2 



(9) 



B. Dimensionless variables 

The problem involves several length scales - particle 
Fermi wavelength X F = tt/pf, the thickness of the layers 
Lx and Li (Lx + Lj = L), the correlation radius of the 
surface inhomogeneities R, and the interface thickness 
D. Another length parameter, the amplitude of inhomo- 
geneities £, is perturbative and enters conductivity as a 
coefficient, 



2e^W_ 

TJ 1 



f(X F ,Li,R,D) 



(10) 



Note, that we consider only the contribution from surface 
roughness and disregard the bulk scattering. As a result, 
the conductivity (|10fl diverges in the limit of vanishing 
inhomogeneities I —* or R — > oo. The proper account 
of bulk scattering [l^ eliminates this divergence. 

It is convenient to measure all the length parameters 
in the units of the Fermi wavelength Xp = n/pp. Instead 
of the interface strength Uo, we use interchangeably two 
equivalent dimensionless parameters g and uq, 



g = uqL/ttXp — 2mUoXpL/Trh 2 



(11) 



(g is convenient for the calculation of the spectrum while 
uo is a proper energy parameter for characterization of 
the conductivity in our setup) . The position of the inter- 
face is characterized by the parameter S, 



5 = L 2 /L. 



(12) 



In computations, S changes from (no second layer) to 1 
(the second layer much wider than the first). It is worth 
repeating that we are looking at the experimental situa- 
tion when the thickness of the first layer is fixed and the 
conductivity is measured as a function of the thickness 
of the second layer (or the overall film thickness). 

The energy spectrum 6j (q) is described by dimension- 
less energy units z%, 



t \ 1 2 



(13) 
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where Zi is given by the solution of the II? Schroedinger 
equation for a quantum well with a 8— type barrier inside 
(see Appendix): 



III. RESULTS 



A. Standard quantum size effect 



sin7rz H — sin (nzS) sin [ttz (1 — 5)] = 0, (14) 



Finally, the conductivity a (L) for the experimental 
setup which has been described above, will be displayed 
by the dimensionless function Jl (L/Xf) , 



2e 2 R 2 



(15) 



for various values of R/Xf, D/Xf, L\/Xf and the 
strength of the barrier u . 

All the figures below present this dimensionless func- 
tion fi,. This function is plotted under the assumption 
that the experiment is performed at fixed thickness of the 
first layer. For uniformity, the figures for weak interfaces 
are plotted for it = 0.1, and for strong interface barriers 
- for uq = 10. The simplest energy spectrum corresponds 
to thin first layers, Xf < Li < 2Xp- Therefore, for trans- 
parency of results, the majority of the data are presented 
for L\/Xf — 1.1 (for comparison, some of the graphs give 
the conductivity for larger Li). 

The computational results below are presented for 
the Gaussian correlation function of the interface 
inhomogeneities, 



C(s) =e 2 exp(-s 2 /2R 2 ) 



(16) 



The angular harmonics for this correlator, which enter 
the transition probabilities in Eq. |JS}, are equal to 



C (0) ( qi , qj )=4irl 2 R 2 [e-Wl (QQ') 
C (1) ( qi , qj ) =4ir£ 2 R 2 \e~ QQ 'h(QQ') 



(17) 



-(Q-Q') 2 /2 



where Q = q^i?, Q' = q jR. 

Analysis of QSE in Ref. for ultrathin films with 
scattering by the film walls demonstrated that the results 
for all types of correlators are qualitatively the same as 
for the Gaussian one when For large inhomo- 

geneities, R 3> L, the results for all types of correla- 
tors with exponential power spectra are similar to those 
for the Gaussian correlator and are qualitatively differ- 
ent from the power-law correlators. The results for the 
power-law correlators are less interesting: such films al- 
ways exhibit the standard saw-like QSE irrespective of 
the value of R because of the wider fluctuations of the 
inhomogeneity sizes. Therefore, in this paper we consider 
only the exponential correlators with a well-defined size 
of inhomogeneities. 



The standard quantum size effect (QSE) in films man- 
ifests itself by a saw-like dependence of the conductivity 
a on the film thickness L [lfl lllj . The positions of the 
singularities - the saw teeth - correspond to the values 
of thickness at which a new energy miniband 6j becomes 
accessible. The amplitude of the conductivity drop in 
such a singular point depends, in the case of scattering 
by surface inhomogeneities, on the effectiveness of the 
roughness-driven interband transitions. If the probabil- 
ity of such transitions Wi^j, is small in comparison to 
the rate of the intraband scattering Wu , the singularities 
in the curves a (L) are almost completely suppressed and 
the standard QSE disappears [9j. 

Analysis of the roughness-driven transition probabili- 
ties for surface scattering in Ref. for different classes of 
surface roughness showed that, when the average size of 
inhomogeneities, R, is much smaller than the film thick- 
ness L, the values of the interband transition probabili- 
ties Wi^j are comparable to that for the intraband scat- 
tering Wa , all the scattering channels are equally impor- 
tant. In this case, the curves a (L) always exhibit the 
standard QSE. The same should be true for scattering 
by the interlayer interfaces. This is illustrated in Figures 
1 and 2 which show a (L) for a weak and strong interface 
potentials uq = 0.1 and uq = 10 respectively. In Figure 1 
the thickness of the first layer is Li = 2.1 Xf, in Figure 2 - 
Li = 1.1 Xp- In both figures, the size of inhomogeneities 
is R = Xp. Both figures exhibit a well- pronounced saw- 
like structure. The positions of the singularities for the 
weak interface are almost equidistant reflecting the fact 
the energy structure is close to that for a square well 
without perturbation inside. The strong interface affects 
the energy spectrum and, therefore, the positions and 
the shapes of the saw teeth. However, at very large film 
thickness L > Li the interface is located very close to 
the well wall and the spectrum start to recover its un- 
perturbed structure. This manifests itself in a recovery 
of the equidistant distribution of the singularities in Fig- 
ure 2 at large L. Because of a peculiar dependence of 
the transition probabilities on the interface strength (see 
Appendix), the conductivity grows much faster with in- 
creasing film thickness in the case of the weak interface 
than for the strong interface. 



B. Quantum size effect for large-scale 
inhomogeneities 

The standard QSE of the type described in the pre- 
vious subsection disappears in the single-layer film when 
the correlation size of inhomogeneities, R, is larger than 
the film thickness, R 3> L, and the correlation function 
in the momentum space £ (q) (the so-called power spec- 
trum of inhomogeneities) decays exponentially at large 
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FIG. 1: Dimensionless conductivity of quantized films, Eq. 
(15), as a function of the film thickness L. The saw-like de- 
pendence is typical for the standard quantum size effect. The 
correlation radius of inhomogeneities R/Xf = 1, the thick- 
ness of the first layer Li/Xp = 2.1, the width of the interface 
d — D/Xf = 0.01, and the strength of the barrier uq — 0.1. 




FIG. 2: Standard QSE in conductivity of quantized films, Eq. 
(15), as a function of the film thickness L for strong interface 
potential, Uo = 10. The correlation radius of inhomogeneities 
R/\f = 1, the thickness of the first layer L\/Xf = 1.1, the 
width of the interface d = 10 -4 . 



wave numbers q. Instead, the single-layer films exhibit 
an anomalous QSE 0. 

The explanation involves the interband transitions. It 
seems that at large R the off-diagonal Wi^k are small and 
the interband transitions are suppressed. However, at 
certain values of large L, few of the elements Wi^k, which 
are close to the main diagonal, could become comparable 



to Wu even for large R. Then the transitions i i + 1 
could become noticeable leading to a drop in conductiv- 
ity. A simple estimate of the peak positions is the follow- 
ing. Scattering by surface inhomogeneities changes the 
tangential momentum by Aq ~ ir/R. This is sufficient 
for the interband transition when this Aq ~ qi — qi+\. 
When the number of occupied minibands is large, the 
lateral Fermi momentum for the gliding electrons, i.e., 
electrons from the miniband with a relatively small index 
i, qi ~ pp. For such electrons, q 2 — q 2 +1 ~ 2nAq/Xp ~ 
2tt 2 /RXf. On the other hand, the energy conservation 
law dictates q 2 — q 2 +1 = (2i + 1) n 2 /L 2 . Accordingly, with 
increasing L the transition channel i «-» i + 1 opens at 
L 2 ~ (i + 1/2) RXf- The opening of a new scattering 
channel in the points 

U ~ y/(i + 1/2) RX F (18) 

is always accompanied by a drop in conductivity. The 
first such drop occurs for the electrons in the lowest mini- 
band ex (q) with i = 1, i.e., for the grazing electrons. 
Note, that these particular electrons contribute the most 
to the conductivity. Since the electrons from the lowest 
miniband are responsible for the dominant contribution 
to the conductivity, the conductivity drops almost by half 
in the point L\ ~ ^SRXf/2 where Wu becomes compa- 
rable to Wu and the effective cross-section doubles. [In 
the quasiclassical film without bulk scattering, the cur- 
rent, which is an integral over momenta, diverges when 
the component of momentum perpendicular to the film 
goes to zero, i.e., for the grazing electrons. Without the 
bulk scattering, the conductivity is finite only because of 
the quantum cut-off at p x = tt/L]. 

The anticipation was that this type of QSE should 
manifest itself also for the interface scattering in mul- 
tilayer films at R 3> L for exponentially decaying surface 
correlators. Indeed, such a picture can be observed in 
Figures 3 and 4 for uq = 0.1; 10 respectively (in both 
figures, L\ — \.\Xf, R — 200Af). The positions of the 
peaks in Figure 3 for the weak interface are close to Eq. 
(|18|l . In the case of the strong interface, the shift of the 
energy levels from those for an "empty" square well is 
much more noticeable and the positions of the peaks in 
Figure 4 deviate from those given by Eq. I|18(l . At large 
values of L, the positions of the peak with strong inter- 
face become close to the points in which the thickness 
of the second layer, L<z = L — L\, rather than the over- 
all thickness L is given by Eq. 1|18|) . The amplitude of 
the anomalous QSE oscillations grows with the increasing 
strength of the interface approaching that for the impen- 
etrable wall. 

Of course, for the inhomogeneities of the intermediate 
size, the picture exhibits the features of both standard 
and anomalous QSE. As it has already been mentioned, 
our numerical examples address the experiment in which 
the size of the inhomogeneities R is fixed while the thick- 
ness of the film L is changing. In general, at the values 
L < R one should see the smooth anomalous QSE oscilla- 
tions with large period, while at L > R one should, on the 
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FIG. 3: Anomalous QSE in conductivity of quantized films, 
Eq. (15), as a function of the film thickness L. The correlation 
radius of inhomogeneities is large R/Xf = 200, the thickness 
of the first layer L\/\f = 2.1, the width of the interface 
d = 0.1, and the strength of the interface barrier ug = 0.1. 




5 10 15 20 25 
L/A F 

FIG. 4: The same as in Figure 3, but for a much stronger 
interface barrier, Uo = 10. 



same curve, see the reappearance of the standard QSE 
with sharper oscillations with period equal to 1 . Roughly, 
the transitions between the regimes occurs when the dis- 
tance between the peaks of the anomalous QSE, Eq. (jl8|l . 
decreases to the value (L;+i — Li) /Xp ~ 1. In principle, 
the reappearance of the standard QSE should be seen 
in Figures 3 and 4 when the computations are extended 
to sufficiently large L. However, the amplitude of the 
standard QSE oscillations on these curves is very small 
and the reappearance of the oscillations is barely notice- 




FIG. 5: QSE in conductivity of quantized films, Eq. (15), as a 
function of the film thickness L for the intermediate values of 
the size of inhomogeneities, R/Xf ~ 3. The thickness of the 
first layer Li/Xf = 1.1, the width of the interface d = O.f, 
and the interface barrier uq — 0.1. At small L, the curve 
exhibits the smooth oscillations of the anomalous QSE with 
a large period, while QSE for large L recovers the standard 
saw-like shape with the period equal to 1. 

able on the scale of the curve. It is much more illustra- 
tive to demonstrate the effect at intermediate values of 
R when both anomalous and standard QSE oscillations 
have comparable amplitude. This is shown in Figure 5 
for R/Xf = 3 and weak interface uq = 0.1. On the left 
side of the graph one can clearly see smooth " new" oscil- 
lations with a relatively large period, while on the right 
side the oscillations recover the sharp saw-like structure 
with the period equal to 1. 

C. Geometric (fractional) quantum size effect 

To exhibit the QSE oscillations of the previous sub- 
section, Figures 3 and 4 were plotted not for the exact 
S— type interfaces (J2J , © but for a somewhat smeared 
(less sharp) interface 

SU = -UqZ (y, z) [5' [x - Lx) + DS" (x - . (19) 

The interface width D can have two origins. If its origin 
is corrugation-related, then the interface width is given 
by the next term of expansion of the interface barrier in 
£ and is characterized by the same parameters £ and R, 
D 2 - (C 2 ). In this case, depending on the correlation 
function, D ~ t or D ~ £ 2 /R. On the other hand, D 
can originate from some " internal" smearing of the inter- 
face and can exist even without surface inhomogeneities. 
In this case, D is a new independent small parameter. 
Note, that here we are interested in the "smearing" of 
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FIG. 6: Geometric QSE in conductivity of multilayer films. 
The same parameters as in Figure 3, except for a much 
sharper interface, d — 10~ 4 . The fractions near the spikes 
give the values of the resonance positions of the interface 
8 = L 2 /L. 

the interface and not in its " fixed" width so that the av- 
erage of the square of the matrix elements of SU over the 
interface starts from D 2 . In Figures 3,4, the interface 
thickness was chosen as d = D/ Xp = 0.1. 

If the interface is thinner, the character of the curves 
changes dramatically. For example, Figure 6 presents 
the conductivity a (L) exactly for the same values of all 
parameters as in Figure 3 except for the interface thick- 
ness which is now d = D/Xp = 0.0001. The difference 
between the two curves is astonishing. 

The conductivity in Figure 6 exhibits two types of 
spikes. The explanation for first type of spikes is the 
following. The scattering-driven transition probabili- 
ties Wij , Eq. (jnj , contain the factor with the correlation 
function £ (q$ — q' j and the geometric coefficients Gij, 
Eq. (|HJ . For exponentially decaying correlators with large 
R ^S> L, the off-diagonal values of the correlation func- 
tion C (q; — q^J with i ^ j are exponentially small in 
comparison with the diagonal ones, £ (q^ — qj). Then it 
is sufficient to analyze only the diagonal elements of the 
matrix Gij, Eq.© : 

G«=W?(Li)*? (ii)- (2°) 

If, accidentally, the S— type interface is positioned in the 
points in which either \I/ ^ (Li) = or {L\) = 0, then 
the coefficient Gu, and, therefore, the transition proba- 
bility Wa, become zero. This, in turn makes the conduc- 
tivity of electrons in the miniband e,, and , therefore, the 
overall conductivity, almost infinite. The cut-off is de- 
termined by one of three factors: 1) exponentially small 
interband transitions; 2) scattering by other defects such 
as impurities, inhomogeneities of external walls, etc.; 3) 




FIG. 7: Geometric QSE in conductivity of multilayer films. 
The same parameters as in Figure 6, except for a wider first 
layer, L\/\p — 4.0 . The fractions near the spikes give the 
values of the resonance positions of the interface S = L2/L. 



smearing of the interface (|19fl that leads to the averaging 
of Gu , Eq. © , H20(l over a finite interval making it non- 
zero. In this paper, for obvious reasons, we are interested 
in the third option. Note, that in the case of scattering 
by external film walls instead of the interlayer interface, 
the coefficients Gij ~ i 2 j 2 are never equal to zero and 
this type of QSE does not exist. 

The first type of spikes corresponds to ^i (Li) = 0. 
The "resonance" positions of the S— type interface are 
universal and do not depend on the potential strength. 
This is true for all rational points 5 = L2/L. Of course, 
the conductivity of the film becomes infinite for this po- 
sition of the interface only if the corresponding mini- 
band ti is occupied. This means that the integer n in 
the denominator of the corresponding fraction 5 — m/n 
should not exceed the number of the occupied minibands, 
n < S — Int [L/Xf]. Indeed, for points <5 = Lij L = m/n 
there is a number of the wave functions ^ (x) of the 
empty well that have nodes in the points x = L\. Since 
the unperturbed homogeneous potential barrier has a 
5— functional form UqS (x — L\), these wave functions 
^i (x) remain the eigenfunctions of the well with the 
unperturbed barrier UqS (x — L±) inside and retain their 
nodes in the points x = L\. Then the corresponding 
diagonal coefficients Gu are zero making the diagonal 
roughness-driven transition probabilities Wu for particles 
from the miniband equal to zero as well. Since the off- 
diagonal transition probabilities are exponentially small 
in R/L ^> 1, the condition Wu = makes the conduc- 
tivity for particles from the miniband e^, and, therefore, 
the overall conductivity exponentially large in R/L 3> 1. 

The structure of the corresponding resonance spikes 
becomes more and more complicated with an increase in 
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L\ when the structure of the minibands and their occu- 
pancy become more convoluted. The simplest structure 
is observed when L\ is between Xp and 2\p as in Figure 
6. In this case, the observed rational spikes correspond to 
the rational numbers of the form 5 = (n — 1) /n and are 
equidistant with the separation L\/Xp. The first spike 
corresponds to the film with 8 = L%/L = 1/2, the sec- 
ond - to 8 = 2/3, the third - to 6 = 3/4, the fourth - 
to 8 = 4/5, and so on. The odd peaks, with the excep- 
tion of the first one, look wider and consist of bigger and 
smaller sub-peaks. The smaller sub-peaks correspond to 
the geometrical resonance with 8 = (n — 1) /n which is 
described above. The bigger and wider sub-peaks have 
a somewhat different nature and are not universal with 
respect to the barrier strength. These sub-peaks will be 
described later. Note, that the peak 8 = 3/4 is so close to 
the first peak from the other series that these two peaks 
are hardly distinguishable. 

When L\ becomes bigger, the first few geometric res- 
onances can be observed at much narrower second lay- 
ers, well before the point 8 = L2/L = 1/2, while the 
density of the resonances become higher. For example, 
Figure 7 presents the conductivity as a function of thick- 
ness for the film with the same parameters as in Figure 
6 except for the thickness of the first layer which is now 
L\ = 4Af- Though the overall distribution of the peaks 
is now much more complicated, the majority can still be 
understood as the ones generated by the eigenfunctions 
of the empty quantum well with the nodes in the posi- 
tions of the barrier. The complexity of the peak structure 
is explained by the fact that at wider first layer L\ more 
minibands are occupied thus allowing a wider selection of 
the rational numbers that determine the peak positions 
8 = L2/ L = m/n. 

The geometric resonances can coexist with the anoma- 
lous QSE of the previous section if the interface is rela- 
tively strong as in Figure 8 for the same configuration as 
in Figure 6 but with much higher value of uq, uq = 10. 
For weak interfaces, the geometric resonances suppress 
the QSE of the previous subsection which gets restored 
only for bigger values of the interface thickness d. This 
graduate disappearance of the geometric resonances can 
be seen when comparing Figure 6 for uq = 0.1, d = 10~ 4 
with Figure 9 (d = 10~ 2 ) and Figure 10 (d = 10" 1 ). Fig- 
ure 10 presents conductivity for the same configuration 
as Figure 3 but in logarithmic scale. In this scale, one 
can see both the wide QSE oscillations of the previous 
subsection and the only surviving geometric resonance 
at 8 = 1/2. 

Above we explained only the narrower, universal ge- 
ometric resonances at 8 = m/n in Figures 6-9. The 
second, non-universal type of spikes has a similar expla- 
nation - zeroes of Ga. According to Appendix, 

G u = 4*2 *f (L,) 

= -JT A t sin2 {KZiS) [g sin (nzrf) + 2z { cos (ttz.S)] 2 , 

(21) 




5 10 15 20 25 
L/A F 

FIG. 8: Geometric and anomalous QSE in conductivity of 
multilayer films for strong interface potential, uo = 10. The 
rest of the parameters are the same parameters as in Figure 6. 
For easier comparison with Figures 4 and 6, the insert gives 
the same data in logarithmic scale. 




5 10 15 20 25 
L/A F 

FIG. 9: Illustration of the effect of smearing of the interface. 
The same curve as in Figure 6, but for a wider interface, 
d = 0.01. 



where Zi (g, 8) is given by the solution of the ID 
Schrocdinger equation i|14|) for a quantum well with a 
8— type barrier inside. The explicit form of the coeffi- 
cients Ai is not important. The factor sin 2 (wZiS) in Eq. 
(|21|l corresponds to 'J' 2 (Li); its zeroes are responsible for 
the geometric resonances with rational 8 = (n — 1) jn. 
There are no other zeroes of sin 2 (nziS). 

However, Eq. (|21|l also contains the factor in the 
square brackets which corresponds to (Li). The 
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FIG. 10: Illustration of the effect of smearing of the interface. 
The same curve as in Figures 6 and 9, but for an even wider 
interface, d — 0.1. Data as in Figure 3, but in logarithmic 
scale. 



simultaneous solution of Eq. I|14l) for the spectrum, 
gsin (irzS) + 2zcos (irzS) — 0, and equation ^ (Li) = 
yield the following equation for the resonance positions 
of the interface: 



sin [nzi (1 - 26)] = 0, 



which is equivalent to 



z i (g,8)(l-25) = k 



(22) 



(23) 



with integer k. First, there is a universal solution 5=1/2 
at k = 0. In this case, when the interface is exactly in 
the middle of the film, both (Li) and ^ (Li) are zero 
(the former with an even index, the latter - with the odd). 
This explains why the geometric resonance with 6=1/2 
is the most stable one with respect to the smearing of the 
interface. 

The rest of the resonances with k ^ are not uni- 
versal. These resonances explain the earlier unaccounted 
for spikes in Figures 6-8. Since the spectrum zi (g, 5) 
is a complicated function of the interface strength and 
its position, the solution of Eq. (|23fl for k ^ is rather 
complicated. We will give the analytic equation for the 
simplest case of \p < L\ < 2\p when Eq. PS^) is equiv- 
alent to 



L 

x7- 



2z n L\ 
— n + 1 Ap 



(24) 



(the only allowed values of k are k = — n +1). For weak 
interfaces g/z n <C 1 this equation can be rewritten as 



1 



2n 



n-1 



-A„ 



sin (irnS) . 

7TZ„ 



Analysis of Eq. (|25|l shows that several first of such res- 
onance positions of the interface are indeed close to the 
odd rational universal resonances as in Figure 6, 

„ 21 - 1 



(25) 



21 ' 

and separate from the universal resonances with increas- 
ing integer I. The very first resonance at S = 1/2 is, 
as it is explained above, exactly the same as the first 
universal resonance. The reason why these non-universal 
resonances are wider and stronger than the universal ones 
described above is still unclear. 



IV. CONCLUSIONS 

In summary, we analyzed QSE in conductivity of mul- 
tilayer films when the main scattering mechanism is the 
scattering of electrons by random inhomogeneities of the 
interlayer interface. Three different types of QSE are 
predicted. 

The first one is a standard QSE with a typical saw-like 
dependence of the conductivity a on the film thickness L, 
a (L). This effect dominates when the correlation radius 
(size) of the interface inhomogeneities R is much smaller 
than the film thickness, R <C L. This effect should be 
observed for all types of the correlation functions of the 
interface roughness. This effect is easily explained by the 
singularities in the electron density of states related to 
the quantization of motion across the film. 

The second type of QSE is explained not by the 
quantization-driven singularities in the density of states, 
but by the anomalies in the cross-section for scatter- 
ing by interface inhomogeneities. This scattering-driven 
QSE replaces the standard saw-like QSE when the cor- 
relation radius (size) of the interface inhomogeneities is 
large, R ~S> L. This type of QSE manifests itself as 
smooth large-scale oscillations on the dependence a (L) 
and should be observed only when the Fourier image of 
the interface correlation function (the so-called power 
spectrum of inhomogeneities) decays exponentially at 
large momenta. The main difference of this QSE from 
a similar effect in scattering by the film walls [3 is that 
the observation of this effect in multilayer film requires 
certain smearing of the interface. 

The third type of QSE is new and is most unusual. 
This effect manifests itself as a set of very narrow and 
high spikes in a (L) and replaces the scattering-induced 
QSE described above when the interface is narrow. The 
finite cutoff in the spikes can be ensured either by some 
other scattering mechanism or by the smearing of the 
interface. 

The spikes are observed only for certain resonance posi- 
tions of the interface. The number of spikes is determined 
by the relation between the thickness of the layers and the 
Fermi wavelength. The resonance positions of the inter- 
face are described. These positions can be split into two 
general classes. Some of these positions are universal and 
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do not depend on the amplitude of the interface potential 
barrier and correspond to the situations when the ratio 
of the layer widths is given by simple rational fractions. 
The integer in the denominator of such fractions does not 
exceed the number of occupied minibands. The remain- 
ing resonance positions of the interface are non-universal 
and depend on the strength of the interface potential. In 
the case of weak interface, some of these non- universal 
positions are close to the universal ones giving the im- 
pression of a split in the conductivity spikes. 

Too small width of the resonance spikes can impede the 
experimental observation of the geometric resonances. 
The width of the resonance spikes increases and their 
height decreases with increasing smearing of the inter- 
face and the resonance spikes gradually disappear. Note, 
that this disappearance of the resonance spikes is related 
not to the widening, but to the random smearing of the 
interface - the widening of the interface, by itself, results 
just in a shift of the spike positions. The width of the 
universal resonances is equal, by the order of magnitude, 
to the width of the smeared interface D, Eq. (|19JI . or, 
in dimensionless variables, d = D/Xp. The width of 
non-universal resonances is somewhat larger and is less 
sensitive to D; the reason is still unclear. Note, that the 
smearing width D can be much smaller than the phys- 
ical thickness of the interface, which in metals is often 
larger than or of the order of the Fermi wavelength Xp. 
In contrast to this, the smearing parameter d = D/Xf 
can be very small. The most stable spike with respect to 
smearing corresponds to the layers of equal width. The 
wide range of possible values of D, which are determined 
either by the roughness with D of the order of I or £ 2 / R 
or by the "internal" smearing of the interface, makes the 
observation of the geometric resonances possible. 

The spikes in the conductivity occur when the scat- 
tering probabilities for electrons in one of the quantized 
minibands become exactly zero. Since scattering proba- 
bilities for scattering by different interfaces add up, the 
spikes in conductivity of multilayer films with many lay- 
ers can be observed only if the scattering probabilities 
for electrons from one miniband become zero simultane- 
ously for scattering by all the interfaces. This can happen 
only if all the interfaces are located in the universal res- 
onance positions corresponding to the rational fractions 
from the same series. Otherwise, the scattering by inho- 
mogeneities of the "non-resonant interface" will curtail 
the contributions from the resonant ones. This imposes 
a restriction on the number of layers for an observation 
of this type of QSE for a film of fixed overall thickness. 

We analyzed the multilayer films under the condition 
that the disruption in the electron spectrum is caused 
only by the interface potential while electron potential 
deep into the layers is the same for all layers. One can 
imagine a different physical situation when the electron 
potential in different layers differ from each other as in 
Ref. [l^- I n this situation the resonance spikes in con- 
ductivity should be observed when the position of the 
interface coincides with one of the nodes in the wave func- 



tion. It is clear that this occurs at least for certain values 
of the interlayer potential difference AU. 

The calculations in the paper are aimed primarily at 
the experimental setup when the lateral conductivity is 
measured as a function of the film thickness at fixed 
thickness of the first layer (fixed position of the buried 
interface) . The main obstacle for the experimental obser- 
vation of the predicted effect is a rather small width of 
the conductivity spikes and their sensitivity to the posi- 
tion of the interface. On the other hand, this sensitivity 
of QSE to the position of the interface may open the door 
for using this effect for precision control of the interface 
positions in multilayer films. This may be very useful for 
better quality ultrathin films without short-range surface 
inhomogeneities |18|| . Recent experiments with controlled 
ultrathin metal films with buried rough interfaces [l| in- 
dicate that the existing experimental setups are sufficient 
for the observation of the predicted quantum size effect. 

Usually, QSE in conductivity of semiconductor films 
is less pronounced than for the metal films. This is ex- 
plained by the smoother distribution of electrons in non- 
degenerate semiconductors. In the absence of sharp drop 
in the distribution at the Fermi energy, singular features 
in the conductivity, which is an integral over the parti- 
cle distribution, tend to be smeared out. However, the 
universal geometric spikes in conductivity, which are de- 
scribed above, are explained by the zeroes in quantized 
electron wave functions on the interface and have nothing 
to do with the electron distribution. Then these spikes in 
conductivity can be the only striking common feature for 
QSE in multilayer metal and semiconductor films. The 
only obstacle for observation of such spikes in semicon- 
ductors could be a relatively large screening radius which 
may lead to an effective smearing of the interface. 

This work was supported by NSF grant DMR-0077266. 

V. APPENDIX: ENERGY SPECTRUM AND 
MATRIX ELEMENTS 

One-dimensional Schrodinger equation for a square 
well with a ^-functional barrier inside has the form 

t/> (x)" + k 2 ip (x) = u S (x - a) if> (x) , (26) 

where 

fc 2 = 2mE/h 2 , u = 2mU /h 2 . (27) 
The wave functions can be written as 

p2 

V>i = tp (x < Li) = \ —Asinkx, (28) 
V L 

p2 

■02 = ip (x > Li) = y —B sin k (x — L) . 
V 

In dimensionless notations of Sec. II. 2, the equation 
on spectrum acquires the form 

sin (ttz) + - sin (nSz) sin [tt (1 - 8) z] = 0, (29) 
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5 = L-i/L < 1/2, kL = Ttz, g = uqL/tt. 

The normalized coefficients in the wave function 1|28|) 
are equal to 



1 



y/S + (1 - 6) t 2 n + t n sin (z„7r) / z n ir ' 
R — A + 



(30) 



where 



in 



sin TtSzn 



sin [ir (1 - 8) z n ] 



(31) 



The explicit expression for the spectrum 1)29(1 can be 
given in the limiting cases of weak and strong potential 
barriers. If the barrier is weak, g/z -C 1, the spectrum is 



z n = n + A„, A„ ~ — sin 2 (wnS) . 



(32) 



In the opposite case of strong interface g — > oo, the 
spectrum decouples into two independent series of levels 
for each layer: 



z ni = ni/S, z n . 2 = n 2 / (1 - 8) 



(33) 



For large, but finite g, the corrections to the spectrum 
(|33|l can be easily obtained by expansion in z n /g8 or 
Zn/g(l - 8): 



"l 



1 



-=■ 1-- 



n 2 
1-8 



1 

1 - - 



(34) 



An important restriction for Eq. (|34[) is that the energy 
levels in each of the layers, Eq. (|33|l . should not be very 
close to each other. In the case of near degeneracy, the 
two close levels, as usual, repel each other with a resulting 
gap equal to 



irgS (1 - 8) 



(35) 



The above equations should be modified if the interface 
is very close to one of the external walls of the well, i.e., 
if cither <5<Clorl-5<l. 



Note, that if the 5-type barrier is located exactly in the 
node of one the wave functions of the empty well, this 
wave function remains the eigenfunction of the well with 
a barrier inside irrespective of the strength of the barrier. 
This means that the energy levels that correspond to such 
wave functions are not shifted by the presence of the 
barrier. 

The matrix elements of the roughness-related pertur- 
bation @) can be calculated with the help of the above 
functions ijSl ty . P ty : 



(36) 



G nm = [tf m (£i) K (Li) + V m (Lx) <&„ (ZOr 
-A A V 2 

g nm = 3 sin (7:z m 5) sin (-Kz n 5) + z m cos (iiz m 5) sin (irz n S) 

+ z n cos (irz n 5) sin (irz m S) . 
The most important are the diagonal matrix elements 

g nn = sin (nz n S) [g sin (nz n 8) + 2z n cos (nz n S)] (37) 
z n sin (irz n 5) sin [7rz n (1 — 25)] 
sin [nz n (1 - 8)} 

Note that the zeroes of the denominator in Eq. ()37|l are 
cancelled out by the zeroes of A 2 , Eq. (|36fl . When the 
interface has a finite width d, the matrix elements acquire 
the following addition: 



2(tot) 



q 2 ' ^ 

>Jnrn 



d 2 (Ag nm y 



where 



Ag nm = 2z n z m cos (nz n S) cos (nz m S) 

sin (nZnS) sin (Trz m S) + g ■ g„ 



iZy. Z^ 



(38) 



The total matrix element, Gnm\ is never zero. This 
means that the term with d 2 , which originates from the 
smearing of the interface, provides a natural cutoff for 
the conductivity. 

In degenerate metal films, of all the energy minibands 
e n only the minibands with n < Int [L/Xp] are occupied. 
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